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Photoemission (PES) and inverse-photoemission spectra (IPES) for the sulphur-passivated InP(OOl) 
surface are compared with theoretical predictions based on density-functional calculations. As a test 
case for our methods, we also present a corresponding study of the better known Si(lll) surface. 
The reported spectra for InP(001)-S agree well with the calculated ones if the surface is assumed 
to consist of a mixture of two phases, namely, the fully S-covered (2 x 2)-reconstructed structure, 
which contains four S atoms in the surface unit-cell, and a (2 x 2) structure containing two S and 
two P atoms per unit cell. The latter has recently been identified in total-energy calculations as 
well as in core-level spectra of S-passivated InP(OOl) surfaces under annealing. The experimental 
IPES for Si(lll)-(2 x 1) is in excellent agreement with the calculations. The comparison of the 
experimental-PES with our calculations provides additional considerations regarding the nature of 
the sample surface. It is also found that the commonly-used density-of-states approximation to the 
photo- and inverse- photoemission spectra is not valid for these systems. 

PACS numbers: 68.35.Bs,73.20-r,79.60.Dp, 85.40.Ls 



I. INTRODUCTION 

The sulphur-treated indium phospide sur- 
face, InP(001)-S, has recently attracted much attention, 
where the effort towards a microscopic understanding of 
the physics has been spurred by its technological impor- 
tance in surface passivation problems.ETd This surface 
has been studied experimentally using a variety of meth^ 
ods, including X-ray absorption near-edge spectroso 
Raman spectroscopy ,0 core- level spectroscopy (CLS)J 
photoemission (PES) and inverse photoemission (IPESI 
spectroscopy^ low-energy electron diffraction-lLEED),E£l 
and scanning tunneling microscopy (STM)O Theoret- 
ically, it has been examined in detail using total en- 
ergy minimization based on density-functional theory 
(DFT).Eia 

Early studies suggested that the InP(001)-S surface is a 
(lxl) structure that follows the bulk (zinc-blende) pat- 
tern, terminated with the sulphur atoms on the bridge 
sites (see for instance Refs. || and |l0|). The actual situ- 
ation, however, has turned out to be much more com- 
plicated. Thus, LEED data on samples annealed at 
350°C appeared to have a clear (2 x 1) pattern.0 How- 
ever, Raman experiments and theoretical calculations 
indicatedQ that the "as-prepared" low-temperature ma- 
terial arranges into a novel (2 x 2) structure with two 
types of sulphur on the surface, as displayed in Fig. 1. 
This structure has four sulphur atoms on the surface 
unit-cell, belonging to two distinct sublayers about 0.2 
A apart; the bottom sublayer is dimerized, i.e., contains 
a S2 "molecule" , while the top sublayer consists of two S 
atoms in a monomer state. We will refer to this surface 



as the S2SS structure. 

Further theoretical calculations, combined with CLS 
measured as a function of the annealing temperature of 
the surface, showed that other surface structures were al- 
lowed within a given thermal window. The most impor- 
tant of these forms via the exchange of surface-sulphur 
atoms with n bulk-phosphorous atoms which migrate to 
the surface.B It was demonstrated that the most-likely 
structure in the 300-400°C annealing range is a (2 x 2)- 
reconstruction containing essentially one type of sulphur, 
with the surface layer containing equal amounts of S and 
P atoms that form tight S-P bonds (Fig. 2) which are 
slightly tilted so that the P-atoms are on top of the sur- 
face. We will refer to this structure as (SP)2- At much 
higher temperatures, all the surface sulphur gets replaced 
by phosphorous, leading to a P-terminated sulphur-free 
InP(OOl) surface. Although the (SP)2 structure is a 
(2 x 2)-reconstruction, it could easily appear to be (2x1) 
if LEED fails to distinguish between S and P atoms. 
Thus, we are lead to a picture wherein the InP(OOl)- 
S surface is actually a system which could contain a 
mixture of S 2 SS, (SP) 2 , and the P-terminated InP(OOl) 
structure, depending on the kinetics imposed by the an- 
nealing conditions. 

The objective of this paper is to examine the reported 
PES and IPES data for the samples annealed at 350°C by 
comparing them with first-principles DFT calculations of 
PES and IPES for the energetically-likely surface struc- 
tures. It turns out that the best fit to the data is obtained 
by assuming that the surface is in a mixed phase consist- 
ing of S2SS and (SP)2 patches. In order to "benchmark" 
our calculations, we present similar calculations for the 
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(2 x l)-reconstructed Si(lll) surface. The spectra for 
Si(lll)-(2 x 1) available in the [literature are for samples 
prepared by cleavage in vacuuoS-n These spectra also con- 
tain various features which could arise from surface im- 
perfections or even the presence of some admixture of 
the more stable (7 x 7) reconstruction. We have exam- 
ined these issues within our analysis of the reported PES- 
spectra of Si(lll), and they serve as a testing ground of 
the theoretical methods used here. 

The plan of the paper is as follows. We first present, in 
Section ||, the essentials of the theory. In Section [II 
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start by discussing our results for the Si(lll) surface and 
cojapare them with the experimental data of Himpsel et 
al.i2 In particular, we emphasize that that the PES and 
IPES spectra are not well accounted for by the density 
of states (DOS) of occupied or unoccupied levels. How- 
ever, good agreement with theory is obtained when the 
calculation is done using the transition matrix elements 
obtained from the Kohn-Sham energy bands. We then 
proceed to a study of the InP(001)-S case and find that 
it is not possible to explain the experimental PES data 
in terms of a single-phase surface; however, we find good 
agreement if the surface is assumed to consist of a mix- 
ture of S2SS and (SP)2 phases which have been identified 
in annealing experiments. 



II. THEORY 

The geometrical details of the S2SS and (SP)2 struc- 
tures described above were determined within the frame- 
work of DFT total-energy minimization calculations; full 
details can be found in Refs. [?] and ^. Given these struc- 
tures, as in our previous studies, we construct supercells 
containing typically 16 atomic layers and the equivalent 
of 5 layers of vacuum. In the case of the Si(lll) slab, 
the top and bottom 4 layers of the slab have the exper- 
imental geometry of the 2x1 reconstruction, while the 
middle eight layers have the bulk-Si geometry. These su- 
percells are then used to calculate the electronic energy 
states (band structure) of these systems. We use, for 
this purpose, the all-electron, fulUpotential-linear-muffin- 
tin-orbital (FP-LMTO) method;t! in this formalism, the 
wavefunctions ip(r) are of the form: 
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(1) 



where n is the band index and k is the electron wavevec- 
tor; t, L, K, and E are, respectively, the index of the 
atomic species in the unit cell, the angular momen- 
tum, the Hankel function expansion parameter (with 
K 2 =— |e|), and the reference energy E of the Bloch basis 

The transition probability between initial and final 
states, ipi(r) and i/>f(r), induced by the electron-photon 
interaction H int is given by the Fermi golden rule as 
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W if = —\(f\ H int I i)\ 2 5{E } - Ei - Jiu) (2) 

where Tim is the photon energy. For PES, electrons with 
final energy Ef will escape from the sample with kinetic 
energy Eki n = Ef — Ep — <p, where <f> is the workfunc- 
tion. The interaction Hamiltonian is A ■ p, where A is the 
vector potential of the photon field. In many cases, it 
is necessary to replace the external field due to the vec- 
tor potential A of the photons by the total field inclusive 
of induced fields using a Sternheimer-type pracedure, as 
shown for instance by Zangwill and Soven.tZI However, 
in our calculations, we restrict ourselves to the spectrum 
given by the simplest form of the Fermi golden rule. The 
velocity form of the matrix was used in the numerical cal- 
culations. The total photoelectron current per unit solid 
angle in PES is then 



if- 
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The expression for the transition probability for IPES is 
identical to PES; the IPES current is given by 
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where, in this case, Ekin = Ei — Ef — <j>- 

Rather than evaluating explicitly the transition matrix 
elements (/ | A-p\ i), it is often assumed, for simplicity, 
that these are constant in the range of the energies con- 
cerned. In this approximation, the PES/IPES intensities 
become proportional to the DOS of occupied/unoccupied 
levels. Another approximation that is made in this con- 
text is to assume that the final electron state in PES 
(or the initial electron state in IPES) is a single plane- 
wave state. In reality, the final state in PES is probably 
more like a superposition of plane- waves dictated by the 
crystal potential and is best represented by the full, self- 
consistent, Kohn-Sham eigenfunction corresponding to 
the energy of the final state. It should be noted that the 
Kohn-Sham eigenvalues are well known to be a poor ap- 
proximation to the quasi-particle energies, but the eigen- 
functions themselves arc generally a good approximation 
to the quasi-particle wave- functions. Thus, in the present 
study, we treat the PES and IPES calculations as sepa- 
rate problems dealing with separate experiments (as they 
indeed are), and do not attempt to relate them to each 
other, or to the experimental surface bandgaps; rather, 
we determine explicitely the Kohn-Sham eigenstates ap- 
propriate to the occupied/unoccupied levels, and use 
these to evaluate in detail the transition matrix elements 
of Eqs. (3) and (4). In fact, the DOS profile is found to 
yield poor approximations to the measured spectra. This 
fact is of some importance since experimental PES /IPES 
results are often presented as a "mapping" of the "band- 
structure" of a material, and theoretical bandstructures 
are often "corrected" to fit in with the PES-dispersions, 
without deconvoluting the effect of the matrix elements 
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and the non-plane wave character of the final (or initial 
) states in PES (or 1PES). 

In our numerical calculations, the (S-functions appear- 
ing in Eq. ^ and related equations will be replaced by a 
Lorentzian functions of width 0.2 eV. In calculating the 
matrix elements, the wavefunction given in Eq. |l| was 
expanded in plane-waves for ease of computation. An 
energy cutoff of 20 Ry was found to be sufficient for the 
spectral range studied, while test calculations using 30 
Ry were performed in order to establish convergence in 
the expansions. Calculations were carried out for bulk 
Si, Si(lll)-(2xl), the S 2 SS and (SP) 2 surface phases, 
and bulk InP. Since we are dealing with admixtures of 
several phases, i.e., surfaces which may only have short- 
range order, calculations are reported only for the T point 
of the surface Brillouin zone. It was verified that increas- 
ing the number of "bulk-like" layers in the simulation 
supercells helped to augment the bulk-like signal in the 
theoretical spectra, arising from the projected bulk DOS 
at the surface's f point. Indeed, since light penetrates 
many atomic layers, the experimental spectra contain, 
to some extent, contributions from the bulk, as well as 
from resonances due to interactions among the surface 
and bulk excitations. In fact, one effect of coupling a 
surface slab to the bulk is to broaden all the peaks. This 
is similar to the effect of the Lorentzian broadening of 
the delta-functions used here. Also, by considering a suf- 
ficient number of bulk-like layers, we ensure that such 
effects are represented in the theory. Thus we checked 
for the convergence of the density of states in each layer 
as a function of slab thickness by considering systems 
containing up to 16 layers. The surface unit cell was 
always taken to be (2x2); that is, no calculations were 
performed for the P-terminated InP surface free of sul- 
phur resulting from high temperature annealing, which 
has a (2x4) structure, or for the Si(lll)-(7x7) surface 
which is also considered in some of the discussions. 



III. RESULTS AND DISCUSSION 

In order to establish the reliability of our approach, 
we first consider the well-known Si(lll) surface as a test 
case. We compare the results available in the literature 
with our theoretical calculations and return to the prob- 
lem of interest here, viz. the InP(001)-S surface. 

A. PES and IPES for the Si(lll) surface 

The Si(lll)-(2x I) surface was extensively studied in 
the early eighties and its spectroscopy played a pivotal 
role in eliminating the carlsLxpodcls in favour of Pandey's 
7r-bonded chain structure £3'EjEj We use here some of 
the published spectra for the nominal Si(lll)-(2x 1) sur- 
face as a benchmark for our calculations. 



Fig. |3j compares our calculations for the IPES using the 
Kohn-Sham transition matrix elements (solid line), and 
DOS for unoccupied levels flashed line), with the exper- 
imental results of HimpselO (squares) for the Si(lll)- 
(2x1) surface. Here we used a computational supercell 
consisting of a 16-layer slab of Si and 5 layers of vac- 
uum. The absolute position of the theoretical curve was 
displaced to conform with the energy scale reported in 
the experiment. The Fermi energy of the (2x1) surface 
is given, experimentally, to be 0.33 eV above the valence 
band maximum and this has been used in positioning the 
experimental spectrum. It is clear that the IPES calcu- 
lation using the simple Fermi golden rule and the full 
matrix elements is in excellent agreement with experi- 
ment, while the corresponding DOS calculation is a poor 
representation of the experimental spectrum. 

In Fig. f|, we present our calculated PES (dashed 
line) and DOS (dotted line) results for the Si(lll)-(2x 1) 
surface, the calculated bulk PES signal (dot-dashed 
lined)_as well as the experimental results of Himpscl 
et allEl We have not subtracted off the secondary- 
electron background in displaying the experimental spec- 
tra. The Si(lll)-(2x 1) samples used in the experiments 
were cleaved under vacuum and believed to be single- 
phase material. However, scanning tunneling microscopy 
(STM) and other surface analysis techniques that are cur- 
rently available were not available then. 

It is clear from Fig. ^ that the PES calculation re- 
covers the peaks at ~ —1.45 and ~ —0.7 eV seen in 
the experiment. (The DOS calculation also has features 
qualitatively similar to the experiment). However, the 
large, broad intensity seen around ~ —3.0 eV, is not re- 
produced by the PES or the DOS calculation. This -3 
eV broad peak (3eVBP) is less prominent in thcpfixper- 
imental PES curves reported in subsequent work£j The 
bulk spectrum at the f-point shows a peak near -3.5 eV 
and -5.5 eV. It is very well known that perfect-clevage 
surfaces of Si are very difficult to obtain. It has been 
suggested that the 3eVBP is a result of angle-integrated 
photoemission from bulk Si, scattered into the normal di- 
rection by surface imperfections.Ej An additional reason 
for getting part of the angle-integrated spectrum super- 
imposed on the angle-resolved spectrum is the effect of 
diffuse light scattering in the optical detection system p£ 
the particular analyser used in this specific experiment .Ej 

How the angle-integration of the bulk signal by sur- 
face imperfections can lead to a 3eVBP can be under- 
stood within the context of another possible explanation 
of the 3eVBP that we have explored. The Si(lll) 2x1 
structure is metastable with respect to the 7x7 recon- 
struction. There is no evidence suggesting that the (2x 1) 
surface reconstructs to the (7x7) form at room tempera- 
ture. However, the energy dumped on the surface by the 
21 eV photon beam (and the secondary electrons etc., 
generated during the PES) may lead to local annealing 
of the surface structure. Such locally-annealed patches 
may exist in a surface which has already been used for, 
say IPES, and then used for obtaining the PES spec- 
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trum. Hence, we have also explored the possibility that 
the surface may contain (7x7) patches, now known to be 
the ground-state of Si(lll). 

In order to explore this possibility, we show in Fig. |B| 
the experimental angle-resolved PES from Si(lll)-(7x7) 
for small angles in the range 0-10 deg. We need to 
consider only small angles here because the 7x7-phascs 
formed by relaxation of the 2x1 phase would be expected 
to have a set of small orientations which averages to zero 
with reference to the original 2xl-surface. Further, the 
f point, common to both (2x1) and (7x7) surfaces, sam- 
ples only photoelectrons that are normal to the surface. 
The averaged signal from the (7x7) surface for small an- 
gles does show the correct behaviour, giving a 3eVBP. 
When the secondary electron-tail is removed, this broad 
band is seen to be a small, but non-negligible, compo- 
nent of the whole spectrum. Hence, one possible expla- 
nation of the 3eVBP near ~ -3.0 eV would be that (7x7) 
patches are present in the cleaved sample believed to be 
(2x1). 

The angle-integration of the bulk signal by surface im- 
perfections is clearly very analogous to the integration 
of the angle-resolved (7x7) data which contain multiple 
features in the -5 eV to -2.5 eV range and in the same 
manner can lead to a 3eVBP spectral feature. This ex- 
planation of the 3eVBP has the merit of not having to 
invoke the existence of 7x7- reconstructed patches in the 
cleaved Si(lll) sample. 

Nevertheless, al these lead us to emphasizes the impor- 
tance of having data from truly well-characterized, clean 
surfaces when relating theory to experiment. Further, we 
note that even if PES and IPES experiments were car- 
ried out using the same sample, it may still be necessary 
to consider local-annealing of the surface and change of 
structure when going from one probe to another £.3 Based 
on the above discussion and Fig. ||, we are led to conclude 
that the sample in IPES - at least the surface region 
probed by the experiment - was essentially a (2x1) re- 
gion. 

B. PES and IPES for the InP(001)-S surface 

In contrast to Si, the InP surfaces are much less un- 
derstood; this is true in particular of InP(001)-S. Obtain- 
ing good STM images or LEED for the as-prepared sur- 
face has been difficult. The LEED from the as- prepared 
surface is very poor and the photoemission is almost 
isotropic, indicative of poor structural order. In our Ra- 
man study of the problem, we recognized these difficulties 
and examined the signal arising from the phonons associ- 
ated with the splitting of the underlying P-layer (just bfe 
low the In-layer) induced by the surface reconstruction.!!^ 
Unlike the outermost-S layer, this inner-P layer was ex- 
pected to be less influenced by surface contamination. 
PES and IPES experiments have been carried out with 
annealed samples characterized using LEED. However, 



our study of the core-level spectra of these samples un- 
der annealing showed that the thermal treatment intro- 
duces new complications by favouring mixed-component 
phases containing a half-S-half-P-terminated surface.B 
Thus, samples annealed at low temperatures for a short 
time (e.g, ~5-10 min.) could easily be an admixture of 
the InP(001)-S 2 SS and InP(001)-(SP) 2 phases. The for- 
mer is the lowest-energy structure for full S coverage, 
while the latter is the structure favoured under annealing 
when bulk-P atoms migrate to replace the S-atoms which, 
in turn, diffuse into the bulk. When higher tempera- 
tures are used in annealing, the completely P-terminated 
InP(OOl) phase which is free of sulphur also becomes pos- 
sible. 

Fig. |^(a) presents our calculation of the PES for both 
the S2SS (full line) and the (SP) 2 (dashed line) recon- 
structions of the InP(001)-S surface. Also included in 
this figure is the spectrum from bulk InP, projected onto 
the surface T point (dotted line) . The experimental data 
of Mitchell et aln are displayed as squares in Fig. ^(b). 
The PES spectra of each phase (and the bulk) were sepa- 
rately calculated and aligned to have the same d-electron 
peak arising from In atoms (which are not significantly 
affected by surface reconstructions). 

It is clear from Fig. ^ that neither the S2SS spectrum 
nor the (SP) 2 spectrum, alone, can account for the ex- 
perimental data, i.e., the measured spectrum evidently 
contains features from both surface phases. In order to 
show this clearly, we display as the full line in Fig. ^(b), 
a composite spectrum obtained by assuming that the re- 
gion of the surface sampled by the light contains a mix- 
ture of S2SS and (SP)2, in proportions of 30% and 70% 
respectively. The composite spectrum was positioned (by 
a rigid shift of the x-axis) to the experimental spectrum 
and scaled to match the intensity of the most prominant 
peak. Clearly, a reasonable fit to the experimental data 
can be obtained in this way. The experimental spectrum 
will depend on details of sample preparation and anneal- 
ing history. However, an invariant property would be 
that a composite theoretical spectrum constructed from 
the spectra of the three main annealing components, i.e, 
S2SS, (SP) 2 and P-terminated InP(OOl), can always be 
found to match the essential features of a given experi- 
mental spectrum. 

We are therefore led ta,conclude that the samples ex- 
amined by Mitchell et alB using PES actually contain a 
mixture of S2SS and (SP)2 phases. This conclusion is 
in agreement with the observation that the width of the 
features in the experimental PES of InP(001)-S is found 
to be much broader than from cleaved GaAs or InP(llO) 
surfaces. c3 

Fig. ^ presents the results of our IPES calculations for 
a sample containing the same mixture of S2SS and (SP)2 
phases as we have just discussed. It is of course not ev- 
ident that the IPES experiments were carried out using 
the same sample that was used for the PES, i.e., with 
no further aging or surface modification, and therefore 
the surface structure need not be identical. Even with 
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identical samples, the regions of the surface probed by 
the photons in PES, and electrons in IPES need not be 
identical. The agreement between the calculation and 
the experiment, in fact, is very unsatisfactory in the re- 
gion between 2.5 and 5.5 eV, where the theory predicts a 
low intensity. However, the fully P-terminated InP(OOl) 
surface is also a possible product of the annealing pro- 
cess and may be a component of the surface. The PES 
and STM oLihe P-terminated InP(OOl) has been stud- 
ied recently.EjO Experiment shows that this phase has a 
prominant broad absorption band just in the region 2.5- 
5.5 eV (dashed curve; Ref. ^). This hypothesis has the 
additional merit of providing a natural explanation for 
the rather high intensity of the IPES signal in the region 
-1.0 to 0.0 eV: there should be very little intensity in this 
gap region for the surfaces phases containing no P atoms, 
as can be seen from the solid curve. The P-terminated 
InP(001)-(2x4) surface, in contrast, has a smaller gap 
than InP(001)-S, and therefore exhibits a sizable signal 
in the -1.0 to 0.0 eV range. This intensity ledge, together 
with the broad band from 2.5 to 5.0 eV, are therefore con- 
sistent with the hypothesis that the sample contains the 
InP(001)-P-(2x4)-reconstructed phase as well. 

IV. CONCLUSION 

We have analysed the published PES and IPES ex- 
perimental results for the S-passivated InP(OOl) surface 
using theoretical calculations for various surface struc- 
tures which are energetically possible. In this analysis we 
have borne in mind the fact that the InP(001)-S surface 
can change its surface composition and structure under 
annealing, as shown by core-level spectroscopy and by 
total-energy calculations. Our conclusion is that the ex- 
perimental results are consistent with a theoretical spec- 
trum obtained from surfaces containing an admixture of 
several possible phases which could occur under thermal 
annealing. Some of the broad features in the experimen- 
tal PES spectrum of Si(lll) were discussed in terms of 
the effect of surface imperfections or local annealing to 
form 7x7 patches. It is felt that further theoretical and 
experimental study of the InP(001)-S surface would have 
to await the availability of better characterized surfaces. 
However, all the presently-available evidence - Raman, 
CLS, PES, IPES etc.- except LEED, seem to be consis- 
tent with our picture of the as-prepared and moderately- 
annealed surfaces being a 2 x 2-reconstruction, essentially 
S2SS for the former and mostly (SP)2 for the latter. 
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FIG. 1. (a) Top view of the fully-relaxed, fully- S-covered 
InP(OOl) S2SS surface. The atom species are identified in 
the z-direction ([001]) layer sequence shown in (b). The x 
and y directions are [110] and [110], respectively. The atomic 
positions (A) are obtained from total energy minimization. 

FIG. 2. (a) Top view of the the (SP) 2 structure at 
half-sulphur/half-phosphorous coverage, (b) Atomic planes 
in the z-direction. The atomic positions are obtained from 
total energy minimization. 

FIG. 3. IPES from the Si(lll) surface. The spectrum 
(PES) calculated using detailed Kohn-Sham matrix elements 
is shown as the full line, while the dashed line gives the result 
of a DOS calculation. 

FIG. 4. PES from the Si(lll) surface. Experimental re- 
sults for the (2x1) structure is shown as a solid line. The 
theoretical results for the 16-layer 2x1 slab (dashed lines) 
and for the bulk Si(lll) crystal (dot-dashed line) are shown. 



FIG. 5. The solid curve labeled "7x7-aver" in the top 
panel is the average over small-angle (0-10 deg.) contribu- 
tions from experimental angle-resolved PES shown in the bot- 
tom panel. 

FIG. 6. Experimental (squares) and theoretical (lines) 
PES from the InP(001)S surface. The secondary electron 
background has been subtracted from the experimental data, 
(a) The individual PES spectra from the two surface phases 
considered in the text, namely S2SS (full line) and (SP)2 
(dashed line), as well as from bulk InP (dotted line), (b) 
The full line is a composite spectrum for a surface containing 
a mixture of 30% S 2 SS and 70% (SP) 2 . 

FIG. 7. Experimental (squares) and theoretical (full line) 
IPES from the InP(001)S surface. Experimental data for 
the P-terminated InP(001)-(2x4), which is free of sulphur 
and produced in higher-temperature anneals, is also shown 
(dashed curve). 
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